function ejer1()

    
    fprintf('\n integral de  e^x  en  intervalo [0,1]:\n ') 

    S = trapecio_compuesta(0, 1, 100,1)

    fprintf('\n integral de polinomimo de taylor de orden 4 centrado en 0  para e^x:\n ') 
    S = trapecio_compuesta(0, 1, 100,2)

    fprintf('\n integral de polinomimo de taylor de orden 8 centrado en 0  para e^x: \n ') 
    S = trapecio_compuesta(0, 1, 100,3)

    fprintf('\n  integral de polinomimo de interpolador de lagrange para e^x  con s  = 4: \n ') 
    S = trapecio_compuesta(0, 1, 100,4)


    fprintf('\n integral de polinomimo de interpolador de lagrange  para  e^x  con s  = 8: \n ') 
    S = trapecio_compuesta(0, 1, 100,5)
    

endfunction

	
function res = trapecio_compuesta(a, b, N,flag)

    h = (b-a)/N;
    sum = 0;
    n = 0;
    s = 0;

    for i = 1:N-1
        x = a + i*h;
        if  flag == 1
            sum += f(x);
        elseif flag == 2
            n = 4;
            sum += f2(x,n);
        elseif flag == 3	
            n = 8;
            sum += f2(x,n);
        elseif flag == 4
            s = 4;
            sum += f3(x,s);
        elseif flag == 5
            s = 8;
            sum += f3(x,s);
        endif

    endfor	

    if  flag == 1
    	res = (h/2) * (f(a) + 2*sum + f(b));
    elseif flag == 2
       	res = (h/2) * (f2(a,n) + 2*sum + f2(b,n));
    elseif flag == 3
       	res = (h/2) * (f2(a,n) + 2*sum + f2(b,n));
    elseif flag == 4
       	res = (h/2) * (f3(a,s) + 2*sum + f3(b,s));
    elseif flag == 5
        res = (h/2) * (f3(a,s) + 2*sum + f3(b,s));
    endif


endfunction
